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Low temperature carrier transport properties in two-dimensional (2D) semiconductor systems can 
be theoretically well-understood within a mean-field type RPA-Boltzmann theory as being limited 
by scattering from screened Coulomb disorder arising from random quenched charged impurities 
in the environment. In the current work, we derive a number of simple analytical formula, sup¬ 
ported by realistic numerical calculations, for the relevant density, mobility, and temperature range 
where 2D transport should manifest strong intrinsic (i.e., arising purely from electronic effects and 
not from phonon scattering) metallic temperature dependence in different semiconductor materials 
arising entirely from the 2D screening properties, thus providing an explanation for why the strong 
temperature dependence of the 2D resistivity can only be observed in high-quality and low-disorder 
(i.e., high-mobility) 2D samples and also why some high-quality 2D materials (i.e., n-GaAs) manifest 
much weaker metallicity than other materials. We also discuss effects of interaction and disorder 
on the 2D screening properties in this context as well as compare 2D and 3D screening functions to 
comment why such a strong intrinsic temperature dependence arising from screening cannot occur 
in 3D metallic carrier transport. Experimentally verifiable predictions are made about the quantita¬ 
tive magnitude of the maximum possible low-temperature metallicity in 2D systems and the scaling 
behavior of the temperature scale controlling the quantum to classical crossover where the system 
reverses the sign of the temperature derivative of the 2D resistivity at high temperatures. 


I. INTRODUCTION 

The observation of a strong apparent metallic temper¬ 
ature dependence of the 2D electrical resistivity in high- 
quality (i.e., low-disorder) semiconductor systems at low 
carrier densities has become fairly routineii^ during the 
last 20 years ever since the first experimental report of 
such an effective metallic behavior in high-mobility n-Si 
MOSFETs^. Typically, the 2D resistivity p{n,T), where 
n is the 2D carrier density and T is the temperature, 
increases with increasing temperature by a substantial 
fraction in the O.IK - 5K regime at “intermediate” car¬ 
rier densities before phonon effects become operational 
at higher temperatures. At very low density, the sys¬ 
tem becomes a disorder-driven strongly localized insu¬ 
lator with an activated (or variable-range hopping) re¬ 
sistivity whereas at high density, the metallic temper¬ 
ature dependence is suppressed with the resistivity be¬ 
ing essentially temperature-independent (except perhaps 
for weak localization effects at very low temperature^ 
which we ignore in the current work). The 2D metal¬ 
lic temperature dependence being of interest here arises 
from intrinsic electronic effects unrelated to phonon scat¬ 
tering (which produces well-known and well-understood 
temperature dependence in the carrier resistivity of met¬ 
als and semiconductors), and thus the low temperature 
transport being discussed in the current work refers to 
the so-called Bloch-Griineisen regime where phonon scat¬ 
tering is strongly suppressed. 

The low-density (“insulating”) and the high (or inter¬ 
mediate) density (“metallic”) transport regimes are sepa¬ 
rated by a crossover density scale ric (sometimes refereed 
to as a critical density although it is really a crossover 


density scale separating an effective metallic phase for 
n > He from a strongly localized insulating phase for 
n < ric) which depends on the sample “quality”, decreas¬ 
ing (increasing) with decreasing (increasing) amount of 
quenched disorder in the system. This low-temperature 
density-driven crossover behavior across ric in going from 
an effective strongly insulating phase (n < ric) to an ef¬ 
fective metallic phase (n > ric), which is sometimes quite 
sharp, is often referred toi as the two-dimensional metal- 
insulator-transition (2D MIT) - a terminology we will 
use in the current work also although in our picture this 
is not a quantum phase transition at all, but is simply a 
sharp crossover from a strongly-localized insulating phase 
to a weakly-localized metallic phase although the weak 
localization behavior may not manifest itself until the 
temperature is unrealistically lowi^ Although a precise 
experimental characterization of the sample quality (i.e., 
the amount of quenched disorder) is challenging because 
of the unknown nature of the impurity distribution^, 
an approximate characterization is provided by the low- 
temperature sample mobility (/i) at high carrier density 
(sometimes referred to as the “maximum mobility”) with 
higher (lower) sample mobility corresponding to lower 
(higher) critical density. Experiments clearly indicate 
that the critical density ric decreases in a particular ma¬ 
terial system (e.g., Si(lOO)-MOSFETs) with increasing 
sample mobility^i^, thus providing a larger range of car¬ 
rier density (n > ric) where the strong metallic tempera¬ 
ture dependence manifests itself, but this dependence of 
the metallic transport behavior on the sample mobility 
does not directly carry over to a comparison among dif¬ 
ferent materials - for example, the metallic behavior is 
strong (weak) for 2D electrons in Si(100)-MOSFETs (n- 
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GaAs) for /i ^ 2 X 10"^ (2 x 10®) cm^/Vs. Thus, the neces¬ 
sary high mobility for the manifestation of strong metallic 
temperature dependence in the 2D transport properties 
depends strongly on the materials system under consider¬ 
ation although in a given 2D system [e.g., Si (100) MOS- 
FETs], the metallicity is typically enhanced with increas¬ 
ing mobility. Clearly, having a high mobility (low disor¬ 
der) is a necessary, but not a sufficient, condition for the 
manifestation of a strong metallic temperature depen¬ 
dence in the 2D resistivity. Similar to the mobility, the 
density and the temperature range for the manifestation 
of the 2D metallic transport is nonuniversal and strongly 
materials dependent although within the same material 
system, the temperature dependence is stronger (weaker) 
with decreasing (increasing) density as long as n > Uc is 
satisfied. For example, in n-Si(lOO) MOSFET (n-GaAs), 
metallicity is observed for n ^ 10^^ (10^) cm“^ in spite 
of the mobility of the GaAs system being typically two 
orders of magnitude higher! 

The current work is focused on analytical understand¬ 
ing of the various materials parameters which are neces¬ 
sary (and sufficient) for the manifestation of the strong 
2D metallic behavior as reflected in the temperature de¬ 
pendent resistivity of 2D semiconductor carriers. The 
theory developed in this article is based on the highly 
successful mean field model of the metallic tempera¬ 
ture dependence in the 2D resistivity as arising from the 
screened Goulomb disorder in the semiconductor through 
the strong temperature dependence of 2D screening. The 
problem is complex even at the mean field level where 
electron-electron interaction is treated entirely through 
static RPA screening of disorder because the total num¬ 
ber of independent physical parameters is large. In addi¬ 
tion to carrier density (n), temperature (T), and mobility 
(/i) mentioned above, transport in 2D systems depends 
also on carrier effective mass (m), background lattice di¬ 
electric constant (k), valley {g^) and spin {gg) degeneracy 
of the 2D materials, various materials parameters char¬ 
acterizing electron-acoustic phonon scattering in the sys¬ 
tem (phonon velocity, Bloch-Griineisen temperature, de¬ 
formation potential coupling, piezoelectric coupling, etc.) 
determining the phonon scattering contribution to the 
electrical resistivity (which is, by definition, temperature 
dependent and must be negligible in order for the screen¬ 
ing induced temperature dependence to be observable), 
and finally the detailed impurity distribution character¬ 
izing the system disorder (with the maximum mobility 
being the minimal parameter defining the system dis¬ 
order). Given this large a set of relevant independent 
parameters affecting 2D transport properties, it seems at 
first hopeless that anything sensible can be stated an¬ 
alytically about the necessary and sufficient conditions 
for the manifestation of 2D metallicity. We show in the 
current work, however, that a few effective parameters 
actually define the theory reasonably well, providing an 
excellent qualitative picture for when and where one ex¬ 
pects the 2D resistivity to manifest a strong metallic tem¬ 
perature dependence. We also present detailed numerical 


transport results for p{T^ n) in several representative 2D 
systems within the RPA-Boltzmann mean field theory in 
support of our qualitative analytical results. 

The rest of the paper is organized as follows. In section 
II we provide a brief comparative discussion of 2D and 
3D temperature dependent screening properties of elec¬ 
tron liquids within RPA to emphasize the physical origin 
of the strong temperature dependence of 2D resistivity as 
limited by scattering from screened Goulomb disorder. In 
section III we present our main analytical arguments de¬ 
riving the conditions for strong 2D metallicity and em¬ 
phasizing the key role of the dimensionless parameters 
qTF/‘^kp and 'T/'Tf, where qTF^ kp, Tp are respectively 
the 2D Thomas-Fermi screening constant, the 2D Fermi 
wave number [kp oc and the Fermi temperature 

{Tp oc n) in determining 2D metallicity. We also provide 
direct numerical results for p{T, n) to support our analyt¬ 
ical results in Section III. In section IV we theoretically 
consider possible corrections to the 2D screening function 
arising from disorder and electron-electron interaction ef¬ 
fects. We conclude in section V with a summary of our 
results, and discussing open questions and possible future 
directions. 


II. TRANSPORT AND SCREENING 


The density and temperature dependent 2D conduc¬ 
tivity limited by disorder scattering is given within the 
Boltzmann transport theory by 


ne^(r) 

m 


( 1 ) 


where the transport relaxation time, (r) = (r(T, n)), is 
defined by the thermal averaging 



with £k = /2m the noninteraction kinetic energy, 

k = |k| the 2D wave number, and f{sk) is the Fermi 
distribution function. In Eq. m, an integral over the 2D 
wave vector k is implied by the summation. The wave 
vector dependent relaxation time r(fc) is given by the 
Born approximation treatment of disorder scattering^i^ 

^ “ ^')l^ (1 - cos6l)<5(efc - efc/), (3) 


where k, k' are the incident and the scattered car¬ 
rier wave vectors (and 9 the angle between them) with 
the 5-function ensuring energy conservation due to elas¬ 
tic scattering by random quenched charged impurities 
with an effective 2D concentration of rii per unit area. 
The carrier-impurity scattering potential is given by the 
screened Coulomb disorder 'u(q) defined as 


u(q) 


^(g) 

e(g)’ 


( 4 ) 
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where v{q) = 27re^/(Kg) is the 2D Coulomb interaction 
(with K the effective background lattice dielectric con¬ 
stant) and e(g), the carrier dielectric screening function, 
is given within RPA by 


e(g) = l-^i;(g)n(g), (5) 


where n(g) is the finite temperature non-interaction 2D 
polarizability function defined by^ 


n(g) = gs9v 

k 


/(gfc) - /(£|k+q|) 
gfc - g|k-fq| 


( 6 ) 


We will not discuss much the theoretical details for the 
RPA-Boltzmann transport theory for disorder scattering 
as provided in Eqs. O-® above since it has already 
been extensively discussed by us in the literature^. We 
note that the actual quantitative theory takes into ac¬ 
count the realistic quasi-2D nature of the semiconduc¬ 
tor system by incorporating appropriate form factors in 
the Coulomb interaction and the Coulomb disorder us¬ 
ing the realistic quasi-2D confinement wavefunctions of 
the 2D carriers. Also, for 3D systems, Eqs. 0-® ap¬ 
ply equally well except, of course, for the wave vector 
k being 3D and integrals in Eqs. ([2]), ®, and dH]) be¬ 
ing three-dimensional with the 3D Coulomb interaction 
being given by 

To understand the role of screening in determining 
2D transport behavior, it is important to realize that 
the most resistive carrier scattering is the 2kp back- 
scattering (i.e., |k — k'l = 2kp) where an electron on 
the Fermi surface gets scattered backward (with a scat¬ 
tering angle 0 = tt) by disorder. Thus, the dominant 
contribution to the temperature dependence of the resis¬ 
tivity at low temperatures comes from the behavior of 
the screening function n(g) around qK, 2kF. Sometimes 
the 2kF scattering is referred to as the scattering from 
Friedel oscillations^ because the singularity structure of 
the polarizability function aX q = 2kF (the so-called Kohn 
anomalyi^) translates to real space Friedel oscillations^ 
of the screened potential. 

This is true both in 2D and 3D, and we, therefore, show 
in Figs. 1 and 2 respectively the calculated temperature 
dependent screening (or equivalently, the noninteracting 
polarizability) function in 2D and 3D in dimensionless 
units [n(g,T) = n(g, T)/n(0, 0)]. A comparison of the 
two figures (Figs. 1 and 2) clearly brings out the key im¬ 
portance of 2kF screening in determining the 2D metallic 
temperature dependence in the disorder-limited carrier 
resistivity, as was already pointed out by Stern quite a 
while ago^. The temperature dependence of the Friedel 
oscillations in the screening clouds around the charged 
impurity centers turns out to be very strong (weak) in 
2D (3D) electron systems as shown in Figs. [1] and [2] here 
and discussed below. 

First, we note that the 2D screening function (Fig. 1) is 
very strongly (going as y/T/Tp) thermally suppressed at 
q m 2kF compared with very weak (going as 


suppression at long wavelength (g = 0). This low- 
temperature thermal suppression of g ~ 2kF screening in 
2D systems is the underlying physical mechanism lead¬ 
ing to the strong metallic resistivity in 2D systemsi^r— . 
We note that the often used long-wavelength screening 
approximation (i.e. the Thomas-Fermi approximation), 
although well-valid at T = 0 since the 2D screening 
function is constant at T = 0 for 0 < q < 2kF hy 
virtue of the constant energy independent 2D density 
of states, fails completely for the calculation of 2D re¬ 
sistivity at finite temperatures since it predicts a very 
weak temperature-dependent 2D resistivity for T <g; Tp 
whereas the full wave vector dep enden t polarizability, 
which includes the anomalous yjT/Tp suppression of 
screening around q ~ 2kp, predicts a strong linear- 
iii-T/Tp increase of the metallic 2D resistivity at low 
temperatures.— This strong temperature-dependence 
of the 2D 2kF screening function is the mechanism under¬ 
lying strong metallicity in 2D semiconductor systems at 
intermediate densities where the value oiT/Tp is not nec¬ 
essarily small leading therefore to a substantial screen¬ 
ing dependent thermal effect. Physically, with increas¬ 
ing temperature, the screened Coulomb disorder, partic¬ 
ularly for the important scattering wavenumbers around 
2kF, is being enhanced strongly due to thermally sup¬ 
pressed screening, leading to an enhanced resistivity due 
to impurity scattering. 

Second, the 3D screening function in Fig. 2 has qualita¬ 
tively different temperature dependence compared with 
the 2D screening function in Fig. 1. In fact, the temper¬ 
ature dependence of the 3D screening function obeys the 
“expected” Sommerfeld expansion behavior in the sense 
that the low-temperature suppression of screening is a 
weak quadratic correction going as 0{T/Tp)"^. This weak 
quadratic temperature dependence applies both for long- 
wavelength Thomas-Fermi screening (g = 0) as well as for 
2A:i?-screening (g = 2kp) implying weak temperature de¬ 
pendence introduced in the 3D resistivity for T/Tp <C 1 
in sharp contrast to the 2D system where the anoma¬ 
lous 0{yjT/Tp) temperature dependence of screening at 
g = 2kF, which violates the Sommerfeld expansion, leads 
to a strong temperature dependence in the carrier resis¬ 
tivity. Thus, the key to understanding the strong metallic 
temperature dependence in the 2D resistivity is the non- 
analytic temperature dependence of the 2D polarizability 
arising from the cusp at 2kp in the non-interacting 2D 
polarizability leading to the failure of the Sommerfeld 
expansion— r— 

For the sake of completeness we quote below the lead¬ 
ing order analytical temperature-dependence of the po¬ 
larizability function in 2D and 3D systems, whereas in 
Figs. 1 and 2 the full numerically calculated polarizabil¬ 
ity is shown for arbitrary temperatures: 


^2D{q = 2.kF,T) = 1 — 


yj(l-y2)C(i) 



(7) 
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FIG. 1. (a) 2D polarizability 'n{q,T) = Ii{q,T)/N]P as a function of wave vector for various temperatures, T = 0, 0.1, 
0.2, 0.5, l.OTp. (b) 2D polarizability as a function of temperature at g = 0. Inset shows II(g = 0, T) at low temperatures. 
The asymptotic form for TjTp <C 1 is given by II(g = 0,T) = [1 — exp[— Tf/T]. (c) 2D polarizability as a function of 
temperature a,t q — 2kF- Inset shows II(g = 2kF,T) at low temperatures. The asymptotic form for T/Tf <S 1 is given by 
n(q = 2kF,T) = 1 — — \/2)^(l/2) Y^r /Tf, where ("(a:) is the Riemann zeta function. 



FIG. 2. (a) 3D polarizability n(g, T) = II(q,r)/A^F^ as a function of wave vector for various temperatures, T — 0, 0.1, 0.2, 

0.5, I.OTf. (b) 3D polarizability as a function of temperature at g = 0. Inset shows II(g = 0, T) at low temperatures. The 
asymptotic form for T/Tf ^ 1 is given by II(g = 0,r) = 1 — ^(T/Tf)^ (red line), (c) 3D polarizability as a function of 
temperature at g = 2kF- Inset shows n(g = 2kF,T) at low temperatures. The asymptotic form for T(Tf <C 1 is given by 
n(g = 2kF,T) = i - -In^] (red line). 


where ('(cc) is the Riemann zeta function. 

n2tt(9 = 0,T) = l-exp(-TF/T). (8) 

n,„(, = o.r) = i-^ (10) 

In Eqs. dB-dini), Tp = Ep/kp is the Fermi tempera¬ 
ture, and Ji 2 D = n(g,T)/iVF^ and HaF) = n(g,T)/iVF^, 
where NjP = = 0,T = 0) and NpP = n 3 pi(q = 

0,T = 0) are the 2D and 3D density of states, respec¬ 
tively. 

Before concluding this section, we emphasize that 
screening is a vital mechanism for 2D semiconductor 
transport because the disorder in the semiconductor envi¬ 
ronment arises primarily from random quenched charged 


impurities whose long-range Coulomb potential must be 
screened for reasonable theoretical results. Thus, within 
a physical mean field approximation, the 2D charge car¬ 
riers (electrons or holes) are scattered from the screened 
Coulomb disorder, and therefore, any strong tempera¬ 
ture dependence in the screening function, particularly 
for 2fcF-scattering which dominates transport at lower 
temperatures, must necessarily be reflected in the 2D re¬ 
sistivity. 

III. THEORY AND NUMERICAL RESULTS 

Having established the importance of 2D screening in 
producing the strong metallic temperature dependence, 
we now analytically derive a number of conditions con¬ 
straining the magnitude of the metallic temperature de¬ 
pendence of 2D transport properties which would ex¬ 
plain the materials dependence of the metallic behavior 
as well as provide reasons for why this metallic behavior 
remained essentially undiscovered (although there were 
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FIG. 3. Schematic p(T) behavior (for various values of 2D 
carrier density n) for low-mobility (a) and high-mobility (b) 
systems. The figure shows the high ric (a) and low ric (b) 
(i.e., rici > nc 2 ) with weak (strong) temperature dependence 
in p{T) in the metallic phase (n > ric) in (a) [(b)] and with 
very similar exponential insulating temperature dependence 
in the localized phase (n < ric). 


occasional hintsi^) until the 1990s in spite of there being 
numerous experimental investigations of 2D semiconduc¬ 
tor transport properties in the 1970s and 1980si^ 

In Fig. [3] we schematically depict the two distinct 
generic experimentally-observed situations for 2D p(T, n) 
with Fig. EKa) and (b) respectively showing the re¬ 
sistivity p{T) for various density (n) in low-mobility 
(high-disorder) and high-mobility (low-disorder) situa¬ 
tions. The only difference between the two situations 
is that one [Fig. Ela)] has a “high” value of Uc (because 
of stronger disorder) whereas the other [Fig. [ll[b)] has a 
“low” value of ric (because of weaker disorder). Thus, 
Figs. EKa) and (b) qualitatively show the respective 2D 
MIT behaviors in the early (< 1995)^ and the present 
(> I995)i days or in low-mobility 2D systemsii and in 
high-mobility systemsi^, respectively. In Fig. Ela) and 
(b) the temperature dependence of p(n, T) is weak and 
strong respectively for n > ric. We mention in this con¬ 
text the seminal importance of the work of Kravchenko 
and collaborators^iiS who first experimentally established 
the connection between the sample quality and the strong 
temperature dependence of the 2D resistivity in the 
metallic (n > ric) phase using low-temperature trans¬ 
port studies in high-mobility (> 10,000 cm^/Vs) Si- 
MOSFETs. Indeed, it is the seminal 1994 — 95 work 
of Kravchenko and collaborators which created the mod¬ 
ern subject of 2D MIT, serving as the temporal milestone 
separating the early days of 2D MIT— [i.e., Fig.EJa)] from 
the present days [i.e.. Fig. EDb)] of 2D MIT—. We em¬ 


phasize that both Figs. 3(a) and (b) manifest essentially 
identical strongly localized insulating phase for n < ric, 
but differ in the temperature dependence of the effec¬ 
tive metallic phase with (older) lower mobility samples 
[Fig. 3(a)] showing little temperature dependence for 
n > ric and the (newer) higher mobility samples manifest¬ 
ing strong metallic temperature dependence [Fig. 3(b)] 
for n > ric- Below we establish that the key to the 
strong metallic temperature dependence of the 2D re¬ 
sistivity (for n > ric) is having (low-disorder-induced) 
low values of the crossover density ric, which makes p{T) 
manifest somewhat complementary temperature depen¬ 
dence {dp/dT > 0 for n > TT-c and dp/dT < 0 for n < ric) 
on two sides of the 2D MIT as depicted in Fig. 3(b). 
On the other hand, for low-mobility samples where ric 
is necessarily high, the metallic phase (for n > ric) does 
not manifest any intrinsic temperature dependence [as 
shown in Fig. 3(a)] except at high enough temperatures 
where phonon scattering effects (ignored in the current 
work) become important. We emphasize, however, that 
at very high (low) density both kinds of samples (low and 
high disorder in Fig. E]) manifest similarly weak (strong) 
temperature dependence. We focus only on the metallic 
(n > ric) phase using the RPA-Boltzmann theory and dis¬ 
cuss the necessary and sufficient conditions for the man¬ 
ifestation of a strong intrinsic (i.e. not phonon-related) 
temperature dependence in the 2D resistivity. The tran¬ 
sition to the insulating phase has been discussed by us 
elsewhereii^ and is not a part of the current work where 
the focus is entirely on the effective metallic regime of 
n > ric- 

To understand how the strong (weak) metallic temper¬ 
ature dependence (for n > ric) correlates with low (high) 
values of ric, we introduce three independent density de¬ 
pendent temperature scales {Tp, Tbg, Tp) which char¬ 
acterize the temperature dependence of the resistivity in 
the metallic phase. These are the electron temperature 
scale defined by the Fermi temperature (Tp), the phonon 
temperature scale defined by the Bloch-Griineisen tem¬ 
perature {Tbg), and the disorder temperature scale de¬ 
fined by the Dingle temperature: 


ksTp = Ep 


fi^kp 

2m 


( Aim \ 

2m \gsgv) 


( 11 ) 


/ 4 ,^^ \ 1/2 

kpTpG = 2tikpVph = 2hvph (- oc (12) 

\9sgvJ 

kBTD=T = ^(—] (13) 

2 \mfiJ 

Here Ep, kp = {Airrijg^gy)^^^, m, Vph, and F are re¬ 
spectively the 2D Fermi energy, 2D Fermi wave vector, 
the carrier effective mass, the phonon velocity, and the 
impurity-scattering induced level broadening (with p as 
the sample mobility). For simplicity, we have defined the 
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level broadening F = h/^T where r is the transport relax¬ 
ation time defining the 2D mobility ^ = a jne = eTim 
with /i being the maximum mobility - in general, the 
broadening F (and therefore the Dingle temperature To) 
is density-dependent through the density dependence of 
mobility which is a complication we ignore for our def¬ 
inition of Td. [We also mention that often the Dingle 
temperature is defined with an additional factor of tt in 
the denominator giving a smaller value for Td in Eq. 
m-] To keep our considerations general, we assume a 
carrier valley degeneracy and a spin degeneracy gs so 
that the total ground state degeneracy is gsgv — gs = ‘2 
in general except in the presence of a strong applied mag¬ 
netic field which could spin-polarize the system making 
( 7 s = 1 whereas = 1 in general except in Si-MOSFETs 
where other values of gy > 1 are possible because of the 
peculiar multi-valley Si bulk conduction band structure. 
The Fermi temperature Tp defines the intrinsic quantum 
temperature scale for the 2D electrons, and when Tp is 
very large (i.e., n very high since Tp oc n), there cannot 
be any temperature dependence in the metallic resistiv¬ 
ity at low temperatures arising from intrinsic electronic 
effects since T/Tp <C 1. Thus, ric needs to be relatively 
low just in order to keep Tp low so that T/Tp is not too 
small for n > ric before phonon effects become significant. 
The Bloch-Griineisen temperature Tbg (oc kp oc ^/n) 
defines the characteristic temperature scale {T > Tbg) 
for phonon scattering effects to become important in the 
2D metallic resistivity. For T < Tbg, phonon effects are 
strongly suppressed, leading to a weak T^-type {p ~ 5—7) 
very high power law in the 2D resistivity arising from 
phonon scattering whereas for T > Tbg, the phonon 
scattering contribution to the 2D resistivity is linear in 
T (which is universally observed in all 2D semiconductor 
systems in the metallic phase for T > 1 — lOiF depend¬ 
ing on the carrier density). Thus, the observation of 2D 
metallic behavior at low temperatures requires T < Tbg 
since trivial phonon scattering contribution to the resis¬ 
tivity for T > Tbg, which is always present, is not the 
issue here. 

This immediately implies that Tp < Tbg is necessary 
for the manifestation of the strong metallic temperature 
dependence in the resistivity since otherwise (i.e., for 
Tp :$> Tbg) the low-temperature (i.e., T < Tp) resis¬ 
tivity will be already dominated by the p ^ T behav¬ 
ior arising from phonon scattering effects applicable for 
T > Tbg- Since Tp oc n and Tbg oc ^/n, the condition 
Tp < Tbg necessitates a low carrier density leading to 
the conclusion that a large ric would lead to the temper¬ 
ature dependence of metallic resistivity (n > ric) being 
dominated by phonon scattering effects. This, in fact, 
typically happened in the 2D systems studied in the 1970s 
and 1980s where phonon effects dominated the metallic 
resistivity [n > ric) suppressing all intrinsic screening- 
induced temperature effects^. Thus, simple dimensional 
considerations of the characteristic electronic (Tp) and 
phononic {Tbg) temperature scales in the problem lead 
to the inevitable conclusion that any strong metallic tem¬ 


perature dependence arising purely from a quantum elec¬ 
tronic mechanism necessitates Tp < Tbg (or ^t least 
Tp ^ Tbg is not allowed), and hence necessarily a low 
ric (so that Tp is not too large even for n > ric). As an 
aside we mention that in 3D metals Tp ~ lO'^K and the 
phonon temperature scale Tbg is replaced by the Debye 
temperature Qd ^ 10^K, so that Tp ^ Qp always. This 
means that the metallic temperature dependence in the 
resistivity arising purely from an electronic mechanism 
is simply impossible in 3D metals at low temperatures 
where phonon effects always dominate down to low tem¬ 
peratures. There can be a weak T^ contribution to the re¬ 
sistivity in 3D metals arising from electron-electron scat¬ 
tering through umklapp processes which cannot happen 
in the 2D semiconductor systems since the umklapp scat¬ 
tering involves very large lattice-scale momentum trans¬ 
fer not of interest in semiconductor transport. 

The role of the disorder-dependent Dingle tempera¬ 
ture Td in the transport problem is rather subtle and 
is relevant at the lowest temperatures T < Tp where 
Tp acts as a lower cutoff suppressing the temperature 
dependence for T <Td. This is because the strong tem¬ 
perature dependence of carrier screening leading to the 
metallic temperature dependence is cutoff for T < Tphy 
impurity disorder effects parametrized by the Dingle tem¬ 
perature. This is because the strong temperature depen¬ 
dence of the 2D polarizability around 2kp is suppressed 
for T > Tp as disorder rounds off the 2kp screeningi^i 
This is discussed in Sec. IV. Thus, Tp explains why the 
metallic temperature dependence for T < Tp < Tbg 
arising from quantum electronic processes does not per¬ 
sist (even in the absence of weak localization which is 
being ignored here) all the way to T = 0 as it would 
for Tp = 0 (and of course, if the electronic temperature 
can be reduced indefinitely which may be an impossibil¬ 
ity). Thus, the temperature dependence of p{T) in the 
metallic phase {n > riy) is bounded from above (by Tbg) 
and from below (by Tp) with the screening induced tem¬ 
perature dependence being strong only in the window 
Td ^ T < Tp < Tbg- Since Tp oc F oc p,~^ where p is 
the characteristic mobility of the system, a large p (i.e., 
low disorder) is necessary in order to keep Tp (as well 
as ric) small so that the temperature dependence of p{T) 
can show up in an appreciable temperature window satis¬ 
fying Tp < T < Tp < Tbg with Tbg > Tp guaranteeing 
that phonon scattering would not play a role in the 2D 
MIT physics. Note that if Tp > Tbg (ke., in highly 
disordered samples) all metallic temperature dependence 
will be totally suppressed. We emphasize that the 2D 
system must be high quality (i.e., low-disorder and high- 
mobility) so that both ric and Tp are small since both ric 
and Tp decrease approximately linearly with increasing 
mobility. One reason that the metallic temperature de¬ 
pendence in the resistivity manifests itself rather strongly 
in Si-based 2D systems even for relatively modest values 
of p (> 10,000 cm^/Vs) is because the effective Tbg is 
rather high in Si because of the high phonon velocity and 
generally weak electron-phonon coupling. 
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The high mobility low-disorder samples of Kravchenko 
et al. (and of others since then in the modern era of 
2D MIT physics) routinely satisfy the constraint Tjj < 
T < Tp < Tbg enabling the observation of the strong 
metallic temperature dependence since Uc and Tp are 
both low in these high-mobility samples whereas the 
older Si-MOSFET samples (before the Kravchenko era), 
where the 2D MIT phenomenon was studied in the early 
days^ i^^’^^ , had high disorder (and low mobility) and con¬ 
sequently high ric (and Tp,) leading to large Tp :$> Tbg 
(as well as large Tp > Tbg) in the metallic phase 
(n > Tic) so that no metallicity could be observed except 
for phonon scattering effects for T > Tbg- Thus, the 
amount of disorder in the sample leading to low or high 
ric (and Tp) is the key to the manifestation of a strong 
metallic temperature dependence in p(T) for n > ric. 

The condition derived above, Tp < Tc < Tbg where 
Tc = Tp (n = Tic), for the manifestation of the strong 
metallicity in the 2D system for n > ric is only a qualita¬ 
tive necessary condition which allows the intrinsic tem¬ 
perature dependence from the 2D screening effect to show 
up in transport properties, but whether such a metal¬ 
lic temperature dependence would actually be a quan¬ 
titatively strong effect or not depends on certain addi¬ 
tional sufficient conditions which we would discuss below. 
These sufficient conditions ensure that the 2A:p’-screening 
is in fact quantitatively significant, not just that it is 
allowed to be present. 

To give a quantitative description underlying the qual¬ 
itative picture discussed above, we borrow (without 
any derivations) from our earlier-obtaine d^d^i^"^ theoret¬ 
ical results providing p{T) in the 2D effective metallic 
phase assuming that the resistive scattering arises from 
screened Coulomb disorder in the system. The quantita¬ 
tive analytical considerations provided below for p(T, n) 
in the metallic phase serve three purposes: (1) They re¬ 
inforce in a concrete manner the qualitative discussion 
given above establishing how the consideration of the 
characteristic temperature scales Tp, Tbg, and Tp (par¬ 
ticularly, their density and mobility dependence) immedi¬ 
ately leads to the conclusion that ric must be small (i.e., 
low disorder and high mobility) for the 2D MIT phe¬ 
nomenon to be associated with a strongly metallic tem¬ 
perature dependence in p{T,n) for n > Uc] (2) they pro¬ 
vide a quantitative understanding of what low (or high) 
ric actually means in a materials-dependent manner, i.e., 
tell us how large can ric be in a specific system (e.g., 
Si-MOSFETs or 2D GaAs systems) and still manifest a 
strongly temperature-dependent p{T) for n> Uc without 
any phonon effects; and (3) they describe how large or 
small Uc should be in going from one 2D system to an¬ 
other (e.g., from 2D Si-MOSFETs to 2D GaAs quantum 
wells) in order for similar metallic temperature depen¬ 
dence to show up in different 2D systems for n > Uc- 

The Boltzmann transport theory gives^ii^iH the fol¬ 
lowing analytical results for the semiclassical Pi{T) in 
2D electron systems at asymptotically low (T <C Tp) 


and high (T ^ Tp) temperatures, respectively 


Pi[T "C Tp) 




1-b 


+ y 


2x T 
1 + xTp 

3/2 


Tf 


O 



(14) 


Pi{T » Tp) 





(15) 

where x = qrp/^kp and y = 2.646[x/(l -I- x)]'^ . In 
Eqs. (03 and (03 . po = p{T = 0) and pi = 
{h/e^)(ni/npx'^) respectively are the impurity-scattering 
induced semiclassical resistivities (hence pi) character¬ 
izing the low and the high temperature limits, and 
Tp = Ep/kp is the Fermi temperature (with n, 
being the respective 2D carrier density and impurity 
density in the system, and kp = (ffiKnj and 
Ttf = gs 9 vFn^ j kT? are the 2D Fermi wave vector and 
Thomas-Fermi screening wave vector, respectively). We 
do not provide the analytical derivations of Eqs. dSD 
and (03 . which can be obtained using Eqs. (P) - (03 as 
shown in Refs, [isl andl^ While Eqs. (HI and m pro¬ 
vide the metallic contributions to p(T) arising from the 
temperature dependence of the screened Goulomb disor¬ 
der, the acoustic phonon scattering by itself contributes 
also to the temperature dependence^i^ given in the high 
(T ^ Tbg) and low (T <C Tbg) temperature limits by 


Pa(T » Tbg) 


Po + Apii 



(16) 


Pa{T <^Tbg) ^ Po + Bph , (17) 

where Tbg = ‘2-hkpVph is the Bloch-Gliineisen tempera¬ 
ture with Vph as the relevant phonon velocity - the con¬ 
stants Aph, Bph depend on the elastic properties of the 
semiconductor—. 

We immediately note that strong metallicity neces¬ 
sitates Tbg > Tp, which means that we must have 
2hkpVph > ffikp/2m, i.e., kp < AmVph/h. Since kp oc 
i/n, the observation of metallicity is a low-density phe¬ 
nomenon restricted to n < Spyrn^Vp^^/nffi = Up where 
phonon effects are suppressed. For n > Up, phonon ef¬ 
fects become relevant for transport. 

In addition, the screening induced metallic tempera¬ 
ture dependence [Eqs. (fT3 and (IT3 ] can only apply for 
n > Uc since, for n < Uy, strong localization induced in¬ 
sulating behavior will dominate (and the metallic theory 
does not apply). Thus, the metallic behavior is only al¬ 
lowed in an intermediate density window ric < n < rip. 
It follows right away that if the sample is so dirty that 
ric > rip, the metallic behavior simply cannot be observed 
in an experimental sample under any circumstance at 
any temperature! Thus, a minimal necessary condition 










for the manifestation of 2D metallic temperature depen¬ 
dence is that 

Uc < Up = Sgym^Vpf^/Trfi^, (18) 

where Uc is the crossover carrier density for the metal- 
to-insulator transition. Since ric obviously increases (de¬ 
creases) with increasing (decreasing) disorder in the sys¬ 
tem, a minimal condition for the observation of metallic- 
ity is that the system must have low disorder or, equiv¬ 
alently, high mobility, at least satisfying Eq. (HI above. 
We also note that Eq. (ITSl) implies [using the Si(lOO)- 
MOSFET materials parameters] an ric < 1-2 x 10^^ 
cm“^ for Si(100)-MOSFETs consistent with experimen¬ 
tal observations in the sense that the modern 2D MIT 
era started with the Kravchenko-Pudalov seminal 1994- 
95 2D transport measurements where the critical den¬ 
sity is indeed less than 10^^ cm“^ whereas the older 
MOSFETs manifested an insulating phase (with acti¬ 
vated conductivity) at a much higher density of ric ^ 
10^^ cm~^ where according to our analysis, no 

temperature-dependent metallic conductivity can be ob¬ 
served except for phonon effects for T > Tbg- 

To see the role of high mobility in the modern 2D MIT 
phenomenon of current interest more clearly we consider 
the specific criterion of the impurity scattering induced 
collisional broadening energy scale defined by the level 
broadening parameter T = ksTo (where Tp is the Dingle 
temperature). Using the loffe-Regel criterion for calcu¬ 
lating Uc^ our condition discussed above, i.e., ric ^ np, 
becomes equivalent to the condition Tp < Tp < Tbg 
for the unambiguous manifestation of the metallic tem¬ 
perature dependence. Using the fact that T = h/2T and 
/i = er/m, we then get the following necessary condition 
on mobility for the manifestation of the metallic phase 

H > he/{2kB'mTp). (19) 

Using a carrier density ric ~ 10^^ cm“^, we get for the 
Si(100)-MOSFETs, g > 21,000 cm^/Vs. For proportion¬ 
ally lower values of ric, the required minimum mobility 
would be proportionally higher, again reinforcing the fact 
that high mobility is a necessary prerequisite for the 2D 
metallic phase (i.e., n > ric) to manifest a strong tem¬ 
perature dependence in the resistivity. It is reassuring 
to note that indeed all modern Si-MOS samples showing 
the canonical 2D MIT behavior after the Kravchenko- 
Pudalov discovery typically have g, > 20, 000 cm^/Vs. 

We note that the above constraints on the critical den¬ 
sity [Eq. dT^ j and the sample mobility (/i) are only the 
necessary conditions, which may not be sufficient for the 
actual manifestation of a strong temperature dependent 
2D resistivity on the metallic (i.e., n > ric) side. For ex¬ 
ample, the actual quantitative screening effect on p(r, n) 
as defined by Eqs. da and m, may simply be too small 
for experimental observation even if the necessary condi¬ 
tion oi ric < n < Up is satisfied. To discuss this is¬ 
sue of sufficient conditions we go back to Eq. a and 
note that for Pi{T) to manifest strong temperature de¬ 
pendence, we must have a; ^ 1 (at least a; > 1) so 


that — ^ ~ 2x/(l -|- a;) is not too small. This requires 
X = qTF/2kp > 1, i.e., qpF > 2kF which translates to 

^ ^ = 2g/^m^e‘^/n'^h'^TT^, ( 20 ) 

where k is the background lattice dielectric constant (as¬ 
suming gs = 2). Obviously n > ric has to be satisfied 
for the 2D system to be in the metallic phase, and so 
metallicity requires the additional sufficient condition of 

nc < riM ^ n, ( 21 ) 

with riM = 2g'/m?‘e^/. For Si(100)-MOSFETs 
with = 2 we get riM ~ 1-2 x 10^^ cm“^, which is 
much larger than ric ~ 10^^ cm“^ for the post-1995 
era 2D MOSFET samples manifesting metallicity in the 
T < Tbg regime of temperatures. This large value of riM, 
however, does explain why low-mobility 2D Si samples do 
not manifest any metallicity since ric is large (> um) in 
such lower quality samples. 

It is gratifying that simple considerations involving Tp, 
Tbg, Tp), and ( 3 'pp/ 2 fcp immediately lead to the predic¬ 
tion that in Si(100)-MOSFETs there would be an ric low 
enough (ric // 10^^ cm“^) for high-mobility (/i > 20, 000 
cm^/Vs) samples to show strong metallic p{T) behav¬ 
ior for n> ric exactly as observed experimentally in the 
post-Kravchenko (> 1995) samples whereas in older low- 
mobility samples with ric 10^^ cm~^, there would be 
no metallic p(T) behavior (except for phonon effects for 
T > Tbg) exactly as seen in lower-mobility MOSFET 
systems^. 

What about other 2D systems such as high-mobility 
2D n-GaAs and p-GaAs systems? Below we briefly dis¬ 
cuss quantitative implications of Eqs. (fHl) - (ITU)) for 2D 
GaAs systems with respect to the 2D MIT phenomena. 

First, 2D n-GaAs has m = O.OTnie, g^ = k, = 13, 
and Vph = 4 X 10^ cm/s in contrast to Si(100)-MOSFETs 
(considered above in depth) which have m = 0.19me, 
gv = 2, K = 12, and Vph = 9 x 10^ cm/s. Applying 
Eqs. (ITTl) - (ITOl) to 2D n-GaAs system, we get 

rip Ri 1.5 X 10^°cm“^; ric < 10^°cm“^; 
p ft! 500, OOOcmVUs; riM « 4 x lO^^cm-^. (22) 

This indicates that one would have to go to very low 
carrier density, way below 10^*^ cm“^, to see any metal¬ 
licity in 2D n-GaAs system. Since Tf{K) Ri 4h in n- 
GaAs where h is the carrier density measured in 10^° 
cm“^, the temperature range (T < Tp) for any possible 
metallic behavior would be well below IK. In addition, 
qTF/'^.kp ~ 0.4/-\/h, which means that qTF/‘2.kF = 1 is 
reached only for n Ri 1.6 x 10® cm“®, implying that ob¬ 
serving strong metallicity (i.e., relatively latge dp/dT) in 
2D n-GaAs would necessitate going to carrier density in 
the range of 1 — 2 x 10® cm“® and T < 100 mK, requiring 
electron mobility of 10^ cm®/Vs. Indeed there is only 
one experimental reporli^^ of observing strong metallic 
behavior in 2D n-GaAs, and it required an ultrahigh mo¬ 
bility of 10^ cm“®/Vs and a sample of very low carrier 
density (~ 10® cm“®) in agreement with our estimates. 
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It is easy to convince oneself using Eqs. (HI - cni) 
and 2D p-GaAs parameters that for GaAs 2D holes, the 
metallic behavior should be routinely observable in sam¬ 
ples with mobilities of 10® — 10® cm^/Vs at carrier densi¬ 
ties around 10^® cm“^. This is indeed the experimental 
situation. 

Thus, we have established in this section why older 
Si-MOSFETs did not see 2D MIT phenomenology: It 
is simply because the sample quality was too low and 
consequently the critical density was too high, making it 
impossible for any screening induced temperature effect 
to manifest itself before the phonon induced tempera¬ 
ture effects show up. It may be worthwhile to obtain 
some rough comparative quantitative estimates for the 
metallic temperature dependence in samples with high 
and low disorder in order to contrast older and newer Si- 
MOSFET samples. We provide such a quantitative com¬ 
parison below for two hypothetical Si-MOSFET samples: 
A (high disorder) and B (low disorder) with high-density 
mobilities of 5,000 cm^/Vs (high ric for sample A) and 
50, 000 cm^/Vs (low ric for sample B) 

Sample A (high disorder) has Uc = 10^^ cm“^, which, 
using Eq. (HI leads to 

(Ap/po)^ < 0.2, (23) 

where Ap = p{Tbg) — p(T = Tu) is the temperature 
induced increase in the metallic resistivity (for n > ric) 
arising from the screening effect. 

Sample B (low disorder) has ric = 10^^ cm“^, which, 
using Eq. HI, leads to 

(Ap/po)b > 1.2. (24) 

Thus, sample A (B) would manifest a less than 20 % 
(more than 120 %) increase in the metallic resistivity (for 
n > ric) between To < T < Tbg arising from screen¬ 
ing effects, clearly establishing that having low (high) 
values of the crossover density ric is the crucial element 
of physics determining strong (weak) metallic tempera¬ 
ture dependence in the system. Since the temperature 
range for metallicity (T < Tbg ^ v^) much smaller 
for the lower-disorder sample B, as it has a much lower 
~ 14Ar compared with ~ 35K in sample A, the 
actual manifested temperature dependence would look 
much stronger in sample B, where p{T) will increase 
by a factor of 2 in the T = 0 — lOAT regime compared 
with only a < 10% increase in p(T) for sample A in the 
same temperature (1 — lOAT) range. This simple esti¬ 
mate shows why older lower mobility MOSFET samples, 
extensively studied in the 1970s and 1980s^ with mobili¬ 
ties around 5, 000 cm^/Vs (or less) never manifested any 
strong metallic behavior because of their relatively high 
values of ric whereas the more recently studied MOS¬ 
FET samples with mobilities above 20,000 cm^/Vs (and 
ric ~ 10^^ cm“^ or less) always manifest strong metal¬ 
lic temperature dependence in its resistivity. The mys¬ 
tery of the so-called strong 2D metallic behavior is thus 
connected directly to the relative magnitude of ric as de¬ 
termined by the 2D sample quality. We do, however. 


mention that some Si(lOO) MOSFET samples with rel¬ 
atively higher mobilities manifested observable metallic 
temperature dependence in the measured resistivity as 
far back as in the early 1980sr^ but this was more of an 
exception since Si MOSEFETs with p > 10, 000 cm^/Vs 
were very rare before 1995. 

Recently, a spectacularly strong metallic temperature 
dependent resistivity was observed^ in Si(lll)-based 2D 
electrons with an unprecedented high maximum mobil¬ 
ity of ~ 200, 000 cm^/Vs. This ultra-high-mobility 2D 
Si(lll) electron system has a valley degeneracy of 6 , 
and manifested almost an order of magnitude increase 
in the metallic resistivity (at n ^ 6 x 10 ^^ cm“^) in 
the T = 0.3 — 4 K in contrast to the 2D Si(lOO) MOS- 
FETs which typically manifest at best a factor of 3 
change in the measured resistivity in a similar temper¬ 
ature window^. The strong observed metallicity in this 
high-mobility Si(lll) system arises from its high valley 
degeneracy 5 ^, = 6 , consistent with the bulk 6 -valley min¬ 
ima electronic structure of Si conduction band leading to 
gy = 2 ( 6 ) in Si(lOO) [(111)] 2D systems. In the context 
of this experimental development^^ it may be worthwhile 
to compare Si(lOO) and Si(lll) 2D systems with respect 
to the various parameters (ric, tibg, ?t-m, x = qTp/'^-kF, 
Tp, Tbg, Tp) defining 2D metallic properties. 

Using Eqs. (ll)-(24) incorporating the materials pa¬ 
rameters (m, K, gy, etc.) for Si(lOO) and Si(lll) 2D 
systems we find n^G ^~ 7.5, ~ 

70, 10, ~ 

~ and /T^p°°'' = 1/3. 

The above considerations show that to obtain the same 
value of ric in Si(lOO) and Si(lll) systems necessitates 
that the Si(lll) system has a much larger (at least by 
3 times) mobility whereas the density range (um) upto 
which the metallicity persists is much higher (by a fac¬ 
tor of 70!) in Si(lll) compared with the Si(lOO) sys¬ 
tem. Most importantly, the large valley degeneracy in 
the Si(lll) system implies an effectively large (an order 
of magnitude larger than Si(lOO) system for the same 
2D carrier density) value of a; = qTp/2kp producing a 
very large value of dp/dT in the metallic phase leading 
to a much stronger metallic temperature dependence in 
the resistivity compared with the Si(lOO) system {gy = 2) 
exactly as observed experimentally. Interestingly, our nu¬ 
merical comparison of the Si(lll) 2D system with the 
Si(lOO) 2D system given above suggests, in agreement 
with the experiment^, the intriguing dichotomy that 
while the critical density (and thus the density range 
where strong metallicity is expected) is higher in the for¬ 
mer, the actual temperature dependent fractional change 
in the resistivity is still considerably higher in the Si(lll) 
system even at this higher absolute density compared 
with the Si(lOO) system because of the large valley degen¬ 
eracy operational on the (111) surface! This establishes 
that one cannot compare ric values between two different 
2D materials to conclude about their relative strength of 
metallicity - although Si (111) has relatively higher ric, 
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Si(100)g^=2 m=0.19me Si(lll)g^=2 m=0.3me 



T (K) T (K) 



FIG. 4. Resistivity as a function of temperature for differ¬ 
ent densities (in the unit of The resistivities are 

calculated with the parameters corresponding to (a) Si(lOO), 
(b) Si(lll) with = 2, (c) Si(lll) with = 6, and the 
same impurity configurations (impurities are located at the 
interface of Si and Si02). The vertical red dashed lines indi¬ 
cate the Dingle temperatures corresponding to the mobility 
fj, — 2 X 10®, 5 X lO'^, 2 X 10^, 5 x 10® cm^/Vs (left to right, 
also Td values are shown at top left corner in the figures in¬ 
cluding a value for 10® cm^/Vs, l.e., the highest temperature 
in each figure). The critical density can be calculated from 
Ec = ksTn = r, i.e., = (3s5tj/47r)(e//rfi), which is inde¬ 

pendent on the effective mass for a given mobility. Thus, for 
given mobilities ^ = 10®, 5 x 10®, 2 x 10^, 5 x 10"*, 2 x 10® 
cm^/Vs, we have the critical densities ric = 48.3, 9.7, 2.4, 
0.97, 0.24x10®° cm“° for = 2, and = 144.9, 29.0, 7.2, 
2.9, 0.72x10®° cm~^ for = 6. For different mobility and 
valley degeneracy ric scales as ric oc psPu/ft. (d) Resistivity 
of a GaAs system as a function of temperature for different 
densities (in the unit of 10®°cm“^), which are calculated with 
the parameters corresponding to GaAs with g^ — 1 and the 
interface impurities between GaAs and GaAlAs. The vertical 
red dashed lines indicate the Dingle temperatures correspond¬ 
ing to the mobility fi = 10®^, 10®, 10® cm^/Vs (left to right, 
also Td are shown at top left corner in the figures including 
values for g. = 10'® and 10® cm^/Vs, i.e., lOK and lOOK, re¬ 
spectively). For given mobilities g = 10®, 10®, 10®, 10®, 10® 
cm^/Vs, the critical densities are given by ric = 0.002, 0.024, 
0.242, 2.42, 24.2x10®° cm“^ for = 1, respectively. Val¬ 
ues of Td (ric) indicate the temperature (density) thresholds 
above which the metallicity behavior may manifest itself as 
discussed in the text. 


it still has stronger metallicity compared with Si (100) 
system. 

To establish the qualitative validity of our analytical 
results presented above, we provide in Fig. 4 our detailed 
numerically calculated results for the 2D resistivity in 
2D Si(lOO), Si(lll) (using both = 2 and 6), and n- 
GaAs (g„ = 1) systems directly using Eqs. (P)-® with 
no additional approximations. It is manifestly clear that 


the numerical results establish that the metallicity is the 
strongest (weakest) in Si(lll) (u-GaAs) system exactly 
as our analytical considerations imply. In Fig. 4, we have 
shown by vertical lines various values of the mobility- 
dependent Dingle temperature (Td) which would cut off 
the temperature dependence, explicitly bringing out the 
fact that the low-mobility samples with high disorder 
(and the associated high Tp values) would not mani¬ 
fest any metallic behavior. An alternative statement is 
that low mobility implies high values of metallic density 
in Fig. 4 applying only for n > tic where the metallic 
temperature dependence is weak. This means that the 
observation of any metallic behavior necessitates low val¬ 
ues of Tic where Fig. 4 shows strong metallic temperature 
dependence 


We conclude this section by summarizing our finding 
for the materials dependence of 2D systems manifesting 
strong 2D MIT behavior (i.e., a strong metallic temper¬ 
ature dependence with large dp/dT > 0 for n > ric). 
We find that Tc defined by Tc = Tp{n = ric) must be 
small enough so that Tc < Tbg for phonon effects to 
be negligible at low temperatures. We also need disor¬ 
der to be small enough so that Tp < Tc, and therefore 
Tp) < Tc < Tbg rnust be satisfied as the necessary con¬ 
dition for the manifestation of 2D MIT. The sufficient 
condition is given by qTp/‘2,kp > 1 (or at least, not too 
small) for n > Uc so that dp/dT is not too small. Using 
the known expressions for the relevant variables Tbg, Tp, 
Td, qTF, and kp we conclude that Pv, rn, and Vph should 
be as large as possible [see Eq. Cl], disorder should be 
as small as possible so that g is large [see Eq. Cl], and 
g/,m^ /should be as larger as possible [see Eq. O], 
implying not only large g^ and m, but also small k. This 
immediately leads to the conclusion that high-mobility 
Si(lll) 2D systems will manifest the strongest 2D MIT 
behavior (since ffu = 6 here, and m is large) whereas 2D 
n-GaAs will have the weakest 2D MIT behavior (since 
m = O.OTwe is the smallest here with g^ = 1 ), and high- 
mobility Si(lOO), 2D p-GaAs, and 2D SiGe systems will 
have strong 2D MIT behaviors since m = 0.19 and (?« = 2 
[for Si(100)[, m = 0.4 and = 1 (for 2D p-GaAs) are 
consistent with strong 2D MIT behavior. It is gratifying 
to know that this material-dependence is exactly what 
is manifested experimentally with high-mobility Si(lll) 
2D systems showing2L2^ the strongest 2D metallic be¬ 
havior and 2D n-GaAs showing the weakest 2D metallic 
behavior—. Of course, if the sample mobility is low so 
that the condition Tp < Tc < Tbg is violated, then 
there would be no 2D MIT behavior at all, as happened 
in almost all low-mobility 2D systems prior to 1995. We 
also mention that if the spin degeneracy is lifted (so that 
Ps = 1 instead of 2 ), for example, by the application of 
a strong parallel magnetic field, then 2D MIT behavior 
is suppressed according to the above considerations, and 
as observed experimentally. Our presented numerical re¬ 
sults (Fig. 4) agree with our analytical results. 
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IV. DISORDER AND INTERACTION EFFECTS 
ON 2D POLARIZABILITY 


The mean-field RPA-Boltzmann screening theory ap¬ 
proach to disorder-limited 2D transport used in our anal¬ 
yses so far ignores the effects of disorder and interaction 
on the screening function itself, and incorporates the 
temperature-induced modification (Fig. [T]) of the finite 
wave number 2D polarizability function as the key phys¬ 
ical mechanism controlling the observed intrinsic metallic 
behavior. The two dimensionless parameters controlling 
the metallic temperature dependence of 2D transport are 
T/Tp and qpp/^kp, both of which should be large (or 
at least, not too small) for the manifestation of metallic- 
ity. This immediately leads to the question of how disor¬ 
der and interaction themselves modify the 2D screening 
function and whether the temperature dependence of the 
resistivity arising from the thermal suppression of 2kF- 
screening in the noninteracting 2D polarizability function 
of the clean system is theoretically robust beyond the ze¬ 
roth order mean field RPA-Boltzmann theory used in our 
considerations. This is of course an important, but also 
a very hard, open question whose answer can at best be 
approximate in any attempted theories since the fate of 
an interacting electron system (either in the continuum 
jellium model of an electron liquid as appropriate for our 
system or in the corresponding Mott-Hubbard-Anderson 
model on a lattice) in the presence of disorder is unknown 
as the problem is a true strong-coupling non-perturbative 
problem (with the notorious fermionic sign problem not 
amenable to large scale computer simulations). 

In this section, which should be considered a continua¬ 
tion of the section II of our article, we present some sim¬ 
ple calculations going beyond the RPA theory of screen¬ 
ing including disorder and interaction effects, and argu¬ 
ing that perhaps disorder and interaction, when they are 
not too strong, would not change the picture qualita¬ 
tively, but our claims and findings in this context are 
rather modest and should be taken as very approxi¬ 
mate attempts toward a long-standing unsolved problem. 
There are alternative (and more ambitious) theoretical 
approachesi^i^Sii^ to the problem (of including disorder 
and interaction in the 2D transport theory) in the con¬ 
text of 2D MIT phenomena which are complementary to 
our work (and which also happen to be much less predic¬ 
tive than our theory - the great advantage of our zeroth 
order theory is its simplicity enabling us to make precise 
quantitative predictions as described in section III of this 
article). We mention that the consideration of interaction 
effects (beyond RPA screening) is not just of academic in¬ 
terest here since the physics of the 2D metallic behavior 
(i.e., the manifestation of a strong temperature depen¬ 
dence in the 2D resistivity) is inherently a low-density 
phenomenon (by virtue of the necessity of a; = qTFl2,kp 
being not too small and ric being low so that Tp is not 
too high) as emphasized in the last section. In fact, we 


can rewrite the dimensionless parameter ‘a;’ as: 

X = qrp/^kp = (25) 

where, = me^ jKy/irn is the dimensionless Wigner- 
Seitz radius characterizing the interaction strength in an 
electron liquid and g = gagv is the total ground state de¬ 
generacy. We note that a large (or not too small) value 
of X, as necessary for strong 2D metallicity, implies that 
Tg cannot be too small which then brings into question 
the quantitative validity of the RPA screening approach 
since should be small (r^ < 1) for the quantitative va¬ 
lidity of RPA. Thus, some justification is needed in ignor¬ 
ing interaction effects in a low-density (r^ > 1) electron 
system where l/r^ basically defines the average number 
of electrons participating in a typical screening cloud. 
The main justifications for our mean-field RPA approach 
(other than its simplicity and predictive power) are (1) 
RPA is empirically known to work well for the quantita¬ 
tive description of many interacting Fermi liquid proper¬ 
ties in 3D metals which typically have Vg > and 

(2) perhaps even more importantly, the 2D MIT phe¬ 
nomenon is primarily a “high-temperature” phenomenon 
where T/Tp cannot be too small for the observation of 
metallic behavior, and as such, interaction effects might 
not be too crucial. In particular, T > Tp (h/tau) is nec¬ 
essary for 2D metallicity to manifest itself, making the 
phenomenon essentially a high-temperature phenomenon 
where quantum correlation effects might be small. Nev¬ 
ertheless, it is necessary to investigate both disorder and 
interaction effects on the 2D screening properties going 
beyond RPA which is what we do below. 

In this section, we theoretically consider the 
influence of random charged impurity disorder— 
and exchange-correlation effects arising from electron- 
electron interactions^^ on the 2D electron polarizabil¬ 
ity (or equivalently the screening function). The non¬ 
interacting static 2D screening function in a clean sys¬ 
tem was first calculated by Stern^ within the random 
phase approximation (RPA). In the current work we con¬ 
sider two separate generalizations of the RPA screening 
theory: Inclusion of impurity scattering effects and in¬ 
clusion of electron-electron interaction effects. In both 
cases, we consider theoretical approximation schemes to 
go beyond the simple RPA theory. The full transport 
theory including both disorder and interaction effects on 
screening, however, remains a formidable open challenge 
for the future well beyond the scope of the current work 
although we cite approximate efforts in this direction by 
other groups using alternative (and highly approximate) 
theoretical technique a^^i^°i^^ 

Our reason for considering (separately) both disorder 
and interaction effects on electronic dielectric screening 
is simply the fact that both become important in the 
low carrier density regime, and thus it is important to 
have some approximate estimate of both corrections to 
the basic RPA screening theory. In this context it is im¬ 
portant to emphasize that we have assumed throughout 
that the system remains homogeneous even in the low 
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(a) I - 



(C) n,= [?) 


(d) r- ^ 

FIG. 5. Set of diagrams used in this calculation: (a) the 
self energy corrections, (b) renormalized Green function, (c) 
static polarizability function (n.y) formed from the renormal¬ 
ized Green function, and (d) impurity ladder vertex correc¬ 
tions. 



FIG. 6. The calculated polarizability as a function of wave 
vector for various scattering rate 7 = 0, 0.1, 0.2, 0.5, l.OEp 
(from top to bottom). 


carrier density regime so that the standard ensemble av¬ 
eraged diagrammatic perturbation theory is applicable 
even in the presence of disorder. This may not, however, 
be true in the low carrier density regime in the pres¬ 
ence of long-range Coulomb disorder where linear screen¬ 
ing itself may fail due to the non-perturbative formation 
of charged impurity-induced electron puddles leading to 
an inhomogeneous density landscape which we have dis¬ 
cussed elsewhere^!. In the current work, we assume that 
the system remains homogeneous throughout and con¬ 
sider disorder and interaction effects on the electronic 
polarizability function diagrammatically. 

We Hrst consider disorder effects on the static 2D 
polarizability. Since Coulomb disorder with its long- 
wavelength divergence must necessarily be screened for 
meaningful results, the inclusion of disorder effects 
on screening involves a nonlinear self-consistent theory 
where the screened Coulomb disorder arising from ran¬ 


dom quenched charged impurities in the background both 
determines screening and is determined by it through 
the renormalized Green’s function. Thus, to calculate 
the screening function in the presence of disorder we use 
the renormalized electronic Green’s function due to the 
electron-impurity interaction. The electronic self-energy 
corrections due to the impurity scattering are obtained in 
the non-crossing multiple-scattering approximation [see 
Fig- El a)]. To calculate the self-energy corrections (E) 
the actual electron-impurity Coulomb interaction should 
be used, but the real calculation with long range Coulomb 
interaction is intractable, particularly because of the non¬ 
linear self-consistency requirement. In this calculation we 
approximate the disorder to be of the short-ranged form, 
voS{x - Xo). 

The static polarizability function (11-^) formed from the 
renormalized Green’s function is obtained from the lad¬ 
der vertex Bethe-Salpeter integral equation which can 
be exactly solved for the short-range disorder model [see 
Fig.EIc) and (d)] 


n^{q)=NF 


dio n((7,w) 

27rj 1 - ^n(g,a;)’ 


(26) 


where Np = gsgvXn/2'K is the 2D density of states at 
the Fermi energy, 7 = 2'KniVQNp with an impurity den¬ 
sity (rii) is the disorder scattering strength at the Fermi 
surface, and n(( 7 , w) is given by 


n(g, UJ) = Y. 

k 

^ 2 1 , F{q,uj) + i 

iF{q,uj) F{q,uj) -V 


(27) 


where F{q,uj) = {2/q)^w + /r — S(a;) — g^/4. The 
chemical potential g is calculated self-consistently so 
that, as 7 is changed, the total density n = gsgv f ^f{k) 
is kept constant, where f{k) is the momentum distribu¬ 
tion function in the presence of disorder and given by 
f{k) = ^xaG~^{k,oj) with G being the electron prop¬ 
agator including disorder scattering effects [Fig.El[a) and 
(b)]. We note that the disorder dependent chemical po¬ 
tential decreases approximately linearly with the scatter¬ 
ing strength [i.e., g{'j)/Ep = I — aj/Ep], This behavior 
is quite a contrast to the temperature dependence of the 
chemical potential, in which the chemical potential de¬ 
creases with temperature exponentially at low tempera¬ 
tures [i.e., g{T)/Ep = I — T/Tpexp{—Tp/T)\. In Fig. [ 6 ] 
we show the calculated static polarizability as a func¬ 
tion of wave vector for various scattering rates. Since 
the chemical potential decreases linearly with disorder 
strength the polarizability at g = 0 also decreases lin¬ 
early with disorder strength. Note that the suppression 
of the polarizability due to thermal effects is exponential, 
i.e., (oc exp[—Tp/T]). As shown in Fig. [Slthe sharp cusp 
at q = 2kp is significantly softened by disorder very simi¬ 
lar to the softening by finite temperature effects shown in 
Fig.[T] Thus, the sharp cusp in the electronic polarizabil¬ 
ity function at q = 2kp is rounded by disorder effects in a 
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FIG. 7. The diagrammatic representation of (a) the ze¬ 
roth (Ho) and (b) the first order (Hi) corrections to the 
polarizability of an interacting system. The wiggled lines 
in (b) indicate the bare Coulomb potential, v{q) — jq. 
(c) The calculated first order correction to the polarizability, 
ni(g) = ni(q)/A^F, as a function of wave vector for various 
temperatures, T/Tf = 0.0, 0.1, 0.2, 0.5, 1.0. 


way similar to thermal effects, and therefore, depending 
on whether temperature or disorder is stronger, screen¬ 
ing will be suppressed either by temperature or by dis¬ 
order. This immediately leads to our physical argument 
in section III for why Tp might cut off the metallicity, 
(i.e., for T < Tp), the metallic temperature dependence 
of 2D resistivity will be suppressed by disorder since Tu 
basically is a measure of the disorder strength. This ex¬ 
plains (at least partially) why T > Tu is necessary for the 
manifestation of the 2D metallic behavior or equivalently 
why 2D metallicity is necessarily a high-temperature phe¬ 
nomenon (or equivalently, why disorder or Tu must be 
very small in order for the 2D effective metallic phase to 
manifest itself) as argued in section III heuristically. 

Now we discuss interaction effect on the 2D polariz¬ 
ability also using a perturbative approach. So far we 
have considered the zeroth order polarizability or the 
non-interacting system [i.e., bare bubble of Fig. [71[a)]. 
In section II we discussed the temperature dependence 
of the zeroth order polarizability. We now obtain the fi¬ 
nite temperature polarizability function by going to first 
order in Coulomb interaction in the diagrammatic per¬ 
turbation theory as shown in Fig. EKb). There are three 
diagrams with one Coulomb line in the first-order inter¬ 


action correction to the polarizability. In Fig. [7l[c) we 
show the calculated total result including all three first- 
order diagrams. We note that although each diagram in 
Fig. [7][b) diverges logarithmically, the sum of the three 
first-order diagrams converges which gives us the finite- 
temperature 2D static polarizability up to leading order 
in Coulomb interaction. At zero temperature we have 
ni ((7 = 0) = Npy/^rg/n, where Np = no ((7 = 0) is 
the 2D density of states and Tg = me^ / {kIt?-J pn) is the 
Wigner-Seitz parameter. In the numerical results shown 
in Fig. 7 we use = 1.8 which corresponds to the elec¬ 
tron density of n = 10^^ cm“^ for the n-GaAs system. 
The zero temperature = 0) shows a very sharp 

peak at q = 2kp, which is a direct consequence of the 
2D characteristic scattering arising at the Fermi surface. 
However, as shown in Fig. [7] this sharp peak at q = 2fci7’ is 
significantly softened (suppression of the Kohn anomaly 
as shown in the zeroth-order polarizability in Fig. 1) by 
finite temperature effects. 

The interaction correction to the finite-temperature 
polarizability presented in Fig. 7 implies a much stronger 
thermal suppression of 2fci?-screening than in the cor¬ 
responding non-interacting case shown in Fig. 1 mainly 
due to the strong exchange-induced enhancement of 2kp- 
screening at zero temperature. Such a strong ther¬ 
mal effect would imply a very strong metallic behavior 
in the 2D resistivity which is not observed experimen¬ 
tally. We also mention that there is not much evidence 
for the strong zero-temperature exchange-enhanced 2D 
Kohn anomaly apparent in Fig. 7(b) at T=0. One pos¬ 
sibility is that higher-order interaction corrections would 
cancel out (at least partially) the strong correction in¬ 
dicated by the first order effects. The other possibility, 
which becomes obvious when we compare Figs. 1(a), 6, 
and 7(c), is that in realistic 2D systems, where both dis¬ 
order (Fig. 6) and interaction (Fig. 7) are present, the 
two effects cancel each other out (at least partially), with 
interaction enhancing the thermal effect (Fig. 7) and dis¬ 
order (Fig. 6) suppressing the thermal effect, leading to 
the original RPA screening result shown in Fig. 1 and 
used throughout our transport theory being a reasonable 
approximation. 

Since the collisional level broadening (disorder) smears 
out the cusp at q = 2kp oi the zeroth order polarizabil¬ 
ity, we also expect the suppression of the sharp peak 
(the Kohn anomaly) in the first order correction of the 
screening function. Thus, the sharp peak in the elec¬ 
tronic polarizability function at q = 2kp is rounded by 
both disorder effects and thermal effects. Our general 
finding is that both thermal effects and impurity scat¬ 
tering effects always weaken 2D screening by rounding 
off the singular behavior at 2kp^ whereas interaction by 
itself could enhance screening. 

Much more work will be necessary for a complete un¬ 
derstanding of both disorder and interaction effects on 
2D transport properties, but our approximate results pre¬ 
sented in this section indicate (at least) the possibility of 
a partial cancellation between the two effects, restoring 
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some confidence in our use of RPA screening in under¬ 
standing 2D transport behavior in realistic semiconduc¬ 
tor structures where both disorder and interaction effects 
are undoubtedly present. 


V. DISCUSSION AND CONCLUSION 

In this work, we have investigated the necessary and 
sufhcient conditions for different 2D semiconductor sys¬ 
tems to manifest strong effectively metallic tempera¬ 
ture dependence in their electrical resistivity arising en¬ 
tirely from intrinsic electronic mechanisms in the Bloch- 
Griineisen temperature regime where phonon effects are 
negligible. Using a physically motivated mean-field 
model of RPA-Boltzmann transport theory, where the 
carrier resistivity arises entirely from scattering off ran¬ 
dom quenched charged impurities in the environment, 
we have argued that the strong temperature dependence 
of 2fci;’-screening in 2D systems could by itself produce 
a metallic temperature dependence in qualitative agree¬ 
ment with experimental findings in many 2D systems. 
We have shown that such a temperature-dependent 2kp- 
screening, and therefore the 2D metallic behavior man¬ 
ifested in the strong temperature-dependent resistivity, 
is intrinsic to 2D systems and does not happen in 3D 
metals or semiconductors which manifest the usual Som- 
merfeld thermal behavior. Using the screening theory, we 
have derived a number of simple analytical results, sup¬ 
ported by direct numerical results, putting constraints on 
the temperature, density, and mobility regimes where the 
2D metallic behavior would be most pronounced and pro¬ 
vided a reasonable explanation for why such a 2D metal¬ 
lic behavior was not observed during the numerous stud¬ 
ies of 2D transport in Si-based 2D systems in the 1970s 
and 1980s. Finally, we considered through leading-order 
perturbative approximations how the inclusion of disor¬ 
der and interaction effects (i.e., going beyond the mean 
field RPA theory) could modify 2D screening properties, 
arguing that the two effects, while being strong individu¬ 
ally, oppose each other qualitatively and quantitatively so 
that RPA may not necessarily be a bad approximation for 
understanding transport properties of 2D semiconductor 
systems. 

The motivation for the current work is easy to state. 
Imagine someone has a 2D sample of some material [e.g., 
a Si(lOO) MOSFET] where the only experimental infor¬ 
mation available is the value of maximum mobility in the 
sample at low temperatures (or just the critical density 
for the system to become a strongly localized insulator 
at low temperature). Can we predict the 2D metallic be¬ 
havior in terms of the density and temperature regime 
where a strong temperature-dependent resistivity would 
manifest itself? Our analyses indicate that the answer to 
this question is affirmative. In particular, the necessary 
condition for the 2D metallicity to manifest itself is that 
the inequality Tp < Tc < Tbg must be satisfied where all 
quantities are defined in Sec. Ill and Tc = Tp (n = Uc), 


whereas the sufficient condition is that the value of the 
dimensionless parameter x = qpp/‘i.kp must not be small 
at the density being used for the resistivity measurement. 
This latter condition implies that the experimental den¬ 
sity n (> Uc) for the observation of the metallicity must 
satisfy the inequality nc < n < um where um is de¬ 
fined in Sec. III. Our analytical considerations immedi¬ 
ately lead to the conclusion that the 2D metallicity would 
be the strongest in the 6-valley degenerate Si(lll) 2D 
system, intermediate in Si(lOO) and p-GaAs 2D systems, 
and by far the weakest in 2D n-GaAs system where one 
would have to go down to n ~ 10® cm“® density for 
observing any appreciable metallic temperature depen¬ 
dence in the resistivity necessitating mobilities around 
10^ cm®/Vs. All of these conclusions are consistent with 
our detailed numerical results and experimental results 
from many laboratories. 

Before concluding, we want to discuss two salient fea¬ 
tures (one occurring at low temperatures, T <C Tf, and 
the other at high temperatures, T ~ Tp) of our screening 
theory with concrete and falsifiable predictions. First, 
the theory [see Eq. (HI] predicts a maximum possible 
magnitude of the 2D metallicity as defined by the di¬ 
mensionless parameter dp/dt where p = p{T)/pQ where 
Po = piT = 0) and t = T/Tp. At very low temperature 
(T < Tp), Eq. m gives: 

dp/dt = 2x/{l + x), (28) 

where x = qTp/2,kp. Assuming a very clean system so 
that Td-CT Tp Tbg (and no complications from 
weak localization corrections ignored in our considera¬ 
tion), the above formula predicts that the dimension¬ 
less temperature derivative of a 2D system must vary 
as ^ a; oc at high density {x 1) and become 

just ‘2’ asymptotically at low density (a; ^ I). The clear 
prediction is that all 2D systems would have a maximum 
possible value of two for the dimensionless metallicity 
as defined by the temperature derivative dp/dtl This 
is a concrete prediction which should be experimentally 
checked for different 2D systems by carrying out trans¬ 
port measurements for the best possible samples (highest 
mobility) at the lowest possible temperatures and den¬ 
sities. To the best of our knowledge, this quantitative 
prediction has never been experimentally tested. Any 
experimentally measured deviation from this predicted 
RPA metallicity would indicate the quantitative impor¬ 
tance of interaction or disorder, providing clear directions 
of future theoretical work in the subject. 

The second prediction is the high-temperature 
quantum-classical crossover predicted around T ^ Tp 
in the theory where the transport behavior changes from 
‘metallic-like’ {dp/dT > 0) to ‘insulating-like’ {dp/dT < 
0) at some density-dependent finite temperature T* with 
the asymptotic analytical properties given respectively in 
Eqs. (fT4l) and (fT5l) for T Tp and T Tp. Provided 
that the condition Tbg ^ Tp is satisfied at the partic¬ 
ular density, this temperature-dependent crossover (at a 
constant carrier density) is generically present in all 2D 
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FIG. 8. Numerically calculated T* jTp as a function of the 
dimensionless screening parameter qpp 1‘lkF for the transport 
results presented in Fig. 4 which are calculated with realistic 
quasi-2D Coulomb form factor effect. The figure shows the 
results for five distinct systems for n-GaAs with <?„ = 1, n- 
Si(lOO) with ( 7 „ = 2, n-Si(lll) with = 2, 4, and 6. We 
note that for realistic parameters and carrier densities the 
calculated T* appears always to fall within 0.3 Tf to STp for 
all existing experimental 2D systems. 


MIT experiments where the measured resistivity p(r, n) 
invariably changes its sign at some high crossover tem¬ 
perature T*{n) with T* being defined by dp/cLT = 0 at 
T = T* at a fixed density n. 

The predictions of Eqs. da and (USD respectively are 
that p(T)—pQ ~ T deep in the quantum regime (T <g; Tp) 
and the conductivity a = 1/p goes as ~ T in the clas¬ 
sical regime T > Tp. The high temperature ‘classi¬ 
cal’ linear in temperature behavior of the 2D conduc¬ 
tivity (for T > T*) has indeed been reported in sev¬ 
eral experiments^, but the topic of high-temperature 
(around T ^ T*) transport in 2D systems has not at¬ 
tracted the attention it deserves perhaps because of the 
uncritical focus of the community on whether the density- 
driven transition around n = ric at low temperature is a 
crossover or a true quantum phase transition. The inter¬ 
esting non-monotonic temperature dependence of p{T) 
around T ^ T*{n) rather clearly establishes that the 
2D metallic phase is an effective high-temperature metal 
rather than a new quantum metal, as such, much more 
experimental work should be devoted to the behavior of 
transport around T ^ T*. 

Our numerical results presented in Fig. 4 show that 
the screening theory consistently gives T* ^ Tp whereas 
experimentally the temperature-induced crossover seems 
to occur over a range of temperature T* ~ Tp/2 — 2Tp, 
depending on the 2D material and the carrier density 
involved. It is interesting that T* always falls within a 
factor of 2 of our predicted quantum-classical crossover 
temperature ^ Tp, but more systematic experimental 


work is necessary in pinning down the behavior of T* as 
a function of the dimensionless variable x = qTp/2,kp 
in different materials system, subtracting out possible 
phonon effects which may not be negligible at higher 
temperaturesUnderstanding the detailed quantitative 
behavior of T* in different systems is particularly im¬ 
portant in view of a competing renormalization group 
based two-parameter scaling theor y^°d° for the metallic 
phase which claims, in contrast to our contention, that 
2D MIT is indeed an interaction-driven quantum phase 
transition with the crossover regime T* < Tp playing 
a crucial role in the quantum coherent scaling proper¬ 
ties of the system. There has even been a claim of an 
experimental verification^ of the scaling theory where 
the resistivity data around the crossover temperature T* 
plays a key role in determining the quantum critical prop¬ 
erties in spite of T* ^ Tp (i.e., an essentially classical 
regime) in the experiment! To contrast our screening 
theory with the scaling theory, it is crucial to under¬ 
stand the detailed behavior of T* in different materials. 
We note, however, that the scaling theory, being a the¬ 
ory involving an expansion in l/pv, should be essentially 
exact for the 6-valley Si(lll) 2D system where the ex¬ 
perimentally observed T* ^ 1.5Tp, in agreement with 
our theory, is in apparent disagreement with the scal¬ 
ing theory where the stringent requirement of T* < h/r 
puts T* < 100 mK in the Si(lll) system in sharp con¬ 
trast with experiments^^. The key difference between 
our theory and the scaling theory is that ours is a bal¬ 
listic “high-temperature” theory where T > Tp is nec¬ 
essary for metallicity (and T* ~ Tp necessary for the 
quantum-classical crossover induced maxima in the re¬ 
sistivity) whereas the scaling theory is strictly a theory 
in the diffusive “low-temperature” (T < Tp) regime so 
that T* < Tp < Tp is necessary. In the experimental 
high-mobility 6-valley Si (111) 2D syste m^^^^ , Tp ~ 100 
mK whereas T* ^ Tp ^ 2 K which appears inconsistent 
with the scaling theory (but perfectly consistent with 
our screening theory), but much more work would be 
necessary to settle this important question. The strong 
metallic behavior of the 2D p-GaAs system is also in¬ 
explicable from the scaling theory perspective since it 
has Pv = 1 whereas in our screening theory, the impor¬ 
tant parameter is gyin which is about the same in two- 
valley Si(lOO) system and one-valley p-GaAs system, in¬ 
dicating similar metallic behaviors in these two 2D ma¬ 
terials in the screening theory. We also point out that 
the 2D metallic phase in our screening theory is an ef¬ 
fective ‘high-temperature’ metal in the ballistic regime 
by virtue of the strong screening induced metallicity in 
the TpG > Tp > T > Tp regime with weak localiza¬ 
tion induced logarithmic insulating behavior necessar¬ 
ily showing up in the ‘low-temperature’ diffusive regime 
(T ^ Tp) although such weak-localization effects may 
not always be easy to see experimentally due to electron 
heating effect necessarily present in semiconductor car¬ 
riers at low temperatures. We mention that the very 
low-temperature insulating upturn in the 2D resistivity 
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has been experimentally seen deep in the diffusive regime 
casting doubts on the scaling theory claim of a generic 
quantum phase transition in 2D interacting disordered 
systems 4ii^ 

The calculation of T*(n), the characteristic ‘high’ (~ 
Tf) temperature where the resistivity derivative changes 
its sign from being metallic-like {dp/dT > 0) to being 
insulating-like {dp/dT < 0), is straightforward, if some¬ 
what tedious, in our screening theory as we need to solve 
the integro-differential equation defined by the condition: 

da{n,T)/dT = Q, (29) 

using the integrals given in Eqs. (l)-(6) defining the tem¬ 
perature and density dependent conductivity, a = 1/p. 
First, we note that it follows directly from Eqs. (l)-(6) 
and Eq. (29), neglecting all nonessential complications 
arising from quasi-2D wavefunction and impurity distri¬ 
bution effects, that T*(n) has an approximate scaling 
form for a given material in the strict 2D limit going as: 

T*/Tf = F(x), (30) 

where F{x = qTF/'^kp) is a scaling function of the di¬ 
mensionless screening parameter x and Tf sets the elec¬ 
tronic energy scale in the problem. Of course in real 
systems, the quasi-2D carrier wavefunctions, the actual 
impurity distribution in the sample, and possible pres¬ 
ence of additional scattering mechanisms (not included 
in our Coulomb disorder model) will lead to deviations 
from the perfect scaling form given in Eq. (30), but an 
approximate scaling form is still expected. In Fig. 8 we 
present our directly numerically calculated T*/Tf as a 
function of the dimensionless screening parameter x for 
the transport results presented in Fig. 4 including re¬ 
alistic quasi-2D Coulomb form factor effects. The cal¬ 
culations carried out for five distinct systems in Fig. 8 
(n-GaAs with = 1, n-Si(lOO) with = 2, n-Si(lll) 
with gy = 2, 4, 6) using realistic density values clearly es¬ 
tablish the rather impressive scaling of T*/Tf as a func¬ 
tion of qTF/‘2kF (~ gyvn /Although these sys¬ 
tems encompass very different effective masses [varying 
from 0.07 for GaAs to 0.3 for Si(lll)] and very differ¬ 
ent valley degeneracies, the characteristic temperature 
T*, when scaled with respect to the Fermi temperature 


(Tp ~ n/gym) of the appropriate system, shows a univer¬ 
sal scaling form with respect to qxF/'^kF- This amazing 
scaling property of T* in the screening theory, which has 
not before been pointed out in the literature, can be di¬ 
rectly experimentally tested by measuring T*{n) in var¬ 
ious 2D systems as a function of temperature. We men¬ 
tion that our finding that T*/Tf essentially falls between 
0.5Tf and 2Tf (as can be seen in Fig. 8) for most 2D ex¬ 
perimental systems, and that T*/Tf tends to be larger in 
the Si(lll) gy = 6 system (around 1.5 — 2Tf)^ compared 
with n-GaAs system {T* ~ O.STp)^ with the Si(lOO) 2D 
system being intermediate is consistent with experimen¬ 
tal results. In addition, our finding of T*/Tf increas¬ 
ing with qTF/‘2kF, which implies that T*/Tf should in¬ 
crease (decrease) with deceasing (increasing) carrier den¬ 
sity in the same sample is also consistent with experi¬ 
ments. It is curious that both our screening theory and 
the interaction-based RG theor y^°d° predict T* to be a 
scaling function, albeit for qualitatively different reasons 
with ours being a ‘high-temperature’ ballistic theory for 
T* (^ Tf) and the RG theory being a ‘low-temperature’ 
(T* <^Tf) diffusive theory. 

Finally, we comment on the possibility of screening- 
induced metallic temperature dependence in 2D 
graphene or other related chiral systems where also 
charged impurity scattering is often the dominant resis¬ 
tive scattering mechanism^. It turns out that the chiral 
nature of graphene completely suppresses 2D backscat- 
tering (i.e. essentially no 2A:F-scattering in graphene), 
and therefore, one does not expecl4^ any screening- 
induced metallic temperature dependence in graphene 
(or 3D topological insulator surface layer) resistivity, and 
in fact, very high-quality high-mobility graphene mani¬ 
fests weakly insulating temperature dependence at low 
carrier densities and low temperatures, which has been 
experimentally verified^. 
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